NEW EFFICIENT STEFFENSEN TYPE METHOD FOR 
SOLVING NONLINEAR EQUATIONS 
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p |. Abstract. In the present paper, by approximating the derivatives 

-<^ I in the Kou et al. [I] fourth-order method by central difference quotient, 

j^^ I we obtain new modification of this method free from derivatives. We 

prove the important fact that the method obtained preserve their order 
of convergence, without calculating any derivative. Finally, numerical 

-^ I tests confirm that our method give the better performance as compare 

|7 j I to other well known Steffensen type methods. 
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t^^ ; 1. Introduction 



^ ■ A large number of papers have been written about iterative methods 

O ' foi' the solution of the nonlinear equations. In this paper, we consider 

^^ i the problem of finding a simple root x* of a function f : D O ^ — y 3fJ, 

i.e. f{x*) = and f'{x*) ^ 0. The famous Newton's method for 

finding x* uses the iterative method 



^ 



C^ , Xn-\-l Xf, 



f'iXnY 

starting from some initial value xo- The Newton's method is an impor- 
tant and basic method which converges quadratically in some neigh- 
borhood of simple root x*. However, when the first order derivative of 
the function f{x) is unavailable or is expensive to compute, the New- 
tons method is still restricted in practical applications. In order to 
avoid computing the first order derivative, Steffensen in [2J proposed 
the following derivative-free method. It is well known that the forward- 
difference approximation for f'{x) at x is 

f{x + h)-f{x) 



fix) 



h 
1 
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If the derivative /'(x„) is replaced by the forward-difference approxi- 
mation with h = f'{xn) i.e. 

./. ^ _ f{Xn + fjXn))- fix) 
^ ^"^"^ ^ fiXn) 

the Newton's method becomes 



Xn+1 Xn 



J yXnj 



f{Xn + f{Xn)) - f{Xn)' 

which is the famous Steffensen's method [2]. The Steffensen's method is 
based on forward-difference approximation to derivative. This method 
is a tough competitor of Newton's method. Both the methods are qua- 
dratic convergence, both require two functions evaluation per iteration 
but Steffensen's method is derivative free. The idea of removing deriva- 
tives from the iteration process is very significant. Recently, many high 
order derivative-free methods are built according to the Steffensen's 
method, see [[3], |1], [5], [6], [7], [8], [9], |T0]] and the references therein. 
Jain [Hj presented the third-order method (Jain Method): 

f{Xn)' 
Vn Xji 



fiXn+ fiXn)) - fiXn)' 

j[Xn) /, IN 

''"^' " """ {f{Xn + f{Xn))-f{Xn)}{f{Xn)-f{yn)y ^'^ 

Dehghan et al. [Oi] proposed two third-order Steffensen type method 
(Dehghan Method I): 

2fixnf 

Vn Xn 



f{Xn + f{Xn)) - f{Xn ~ f{Xn)) ' 
_ 2fiXn){fiXn)+fiyn)} ,^ .. 

" f{Xn + f{Xn))-f{Xn-f{Xn)y ^ ' 

and (Dehghan Method II) 

2/(x.)' 



Vn Xn ~r 



f{Xn + f{Xn)) - f{Xn " f{Xn)) ' 



^ ^ _ 2/(x„){/(|/„) - fjXn)} ,^ g. 

fiXn + fiXn))- f{Xn- f{Xn))' 
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Again Dehghan et al. [lO] introduced a new third-order Steffensen 
type method (Dehghan Method III): 

2/(x.)' 



•^n ~r 



f{Xn + f{,Xn)) - f{Xn - f{Xn)) ' 



2/(Xn) 

"^"^^ '^'^ f{yn)f{u) + f{Xr.)f{vy ^'-^^ 

where f{u) = f{xn + /(a:„)) - f{xn - /(a;„)) 
and /(f) = f{yn + /(y„)) - fiVn - fiVn))- 

Recently Cordero et al. [3] presented a fourth-order Steffensen type 
method (Cordero Method): 

2/(x„)2 



Xn 



f{Xn + f{Xn)) - f{Xn - f M) '' 



^^^^ ^ ^^ 2/(^n)' .iMll/M/.1.5) 

f{Xn + f\Xn)) - f{Xn - f (Xn)) ' ^fiy^) - f\Xn) 

Other Steffensen type methods and their applications are discussed in 
[0], 0, |6], H]. 

The purpose of this paper is to develop a new fourth-order derivative- 
free method and give the convergence analysis. This paper is organized 
as follows. In Section 2, we present a new two-step fourth-order itera- 
tive method for solving nonlinear equations. In this method we approx- 
imate the derivative of the function by central difference quotient. The 
new method is free from derivative. We prove that the order of con- 
vergence of the new method is four. Numerical examples show better 
performance of our method in section 4. Section 5 is a short conclusion. 

2. Development of the method and analysis of 

convergence 

Let us consider the fourth-order of method proposed by Kou et al. 
in [I]: 

f{Xn) 



tin Xj-, 



fiXuV 



where 6 & R. The convergence order of these methods is at least 
three, and for 6 = —1 the order is four. Specifically, the third order 
method obtained taking ^ = in (12. ip . and the forth order method 
corresponding to 9 = —1. Here we are going to use this fourth-order 
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method of this family. We try to use forward-difference approximation 
of the derivative appeared in f l2.ip . then the the order of convergence 
of the new method goes down to three and error expression for it is 
(— C2 — C2/ci)e^ + O(e^), where c, is defined later. For this reason, 
we have used the central-difference approximation in (12. ip . obtaining a 
modified Kou et al. method that preserves the convergence order and 
is derivative free. 

By using a symmetric difference quotient 

By approximating the derivative by the central-difference in (12. ip for 
6 = —1, we obtain a new method free from derivatives, that we call 
the modified Kou method free from, derivative (MKDF): 

2/(x„)2 



^n 



f{Xn + f{Xn)) - f{Xn " f{Xn)) ' 

2/(x.){/(x„) + /(!/„)} 



Xn+1 Xn + 

jyXn -h J\,Xn}) - J[Xn - J[Xn>> 

-(2.2) 



f{Xn + f{Xn)) - f{Xn - f M) 
^fiXnf 



{f{Xn + fiXn)) - fiXn - f iXn))}{f (Xn) - f {Vn)} 

Now we are going to prove the method MKDF have orders of conver- 
gence four. 

Theorem 2.1. Let x* E I he a simple zero of a sufficiently differen- 
tiable function f : I 'O^ ^ ^ in an open interval I. If Xq is sufficiently 
close to X* , then the modified Kou method free from derivative defined 
by (12. 2p has order of convergence four and satisfies the error equation 

(ED- 

Proof. By applying the Taylor series expansion theorem and taking 
account f{x*) = 0, we can write 

fixn) = Cien+C2el+Csel+C4e'^+C5el+CGel+C7el+Csel+0{el), (2.3) 

where Ck = [^ , A; = 1, 2, ... and e„ be the error in Xn after n iterations 

I.e. Cji Xfi X , 

f{Xn + f{Xn)) = (Ci + Ci)e„ + (C2C1 + 3C2Ci + €2)6^ 

+ (C3C? + 3c3C? + 2cicl + 4c3Ci + 2cl + C3)el 

+ (c4 + C2{cl + 2C1C3) + 5CiC4 + 5C2C3 + Qcjci 

+4ciC4 + C1C4 + 6C1C2C3 + 3c^C2C3)e^ 
+ 0{el). (2.4) 
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and 

f{Xn - f{Xn)) = (-Ci + Ci)e„ + (C2C1 - 3C2Ci + €2)6^ 

+ (-C3C? + 3c3C? + 2cicl - AC3C1 - 2cl + C3)el 

+ (c4 + C2{cl + 2C1C3) - 5CiC4 - 5C2C3 + Qcjci 

-Ac\c4^ + c\ci + 6C1C2C3 - 3clc2Cs)e'^ 
+ Oiel). (2.5) 

Further more it can be easily find 

2.fM' 

f{Xn + fiXn)) - f{Xn " /{Xn)) 

= en + (-C2/ci)e2 + {{2cl)/cl - (2c3)/ci - CiOi)el 

+ (C2C3 - 4ciC4 - (3c4)/ci - (4c^)/c? + (7C2C3)/C?)e^ 

+0(e^). (2.6) 

By considering this relation and expression of Hn in the equation (12. 2p , 
we obtain 

Vn = x* + {c2/ci)el - {{2cl)/cl - (2C3)/Ci - CiC3)e3 

-(C2C3 - 4ciC4 - (3c4)/ci - (4c^)/c? + (7c2C3)/c?)e^ 

+0{el). (2.7) 

At this time, we should expand /(?/„) around the root by taking into 
consideration (12. 7p . Accordingly, we have 

fiVn) = C2el + Ci{ciC3 + {2c3)/ci-{2cl)/cl)el 

+ (ci(4ciC4 - C2C3 + (3c4)/ci + {AcD/cl - {7c2C3)/cl) + cl/cDef^ 
+Oiel). (2.8) 

By using (Q, UM, (D and ([23]) in the last expression of (Q, we 
obtain 

Cn+l = (-4CiC4 - 3C2C3 - (llC4)/Ci - (2c^)/c!^ - (17 C2C3) / cj) 6^ 

+0(e^). (2.9) 

D 
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3. Numerical Tests 

In this section, in order to compare the our new method with Stef- 
fensen method, Jain method, Dehghan method I, Dehghan method II, 
Dehghan method III and Cordero method, we give some numerical 
examples. For this consider the following functions: 

Table 1. Test functions and their roots. 



Non-linear functions 



Roots 



flix) 


= sin^(a;) - x^ + 1 


1.404492 


f2{x) 


= x^-e^ -3x + 2 


0.257530 


h{x) 


= cos{x) — X 


0.739085 


Mx) 


= cos(a;) — xe^ + x'^ 


0.639154 


Mx) 


= e"^ — 1.5 — arctan(x) 


0.767653 


Mx) 


= 8a; — cos(a;) — 2x'^ 


0.128077 


frix) 


= ln(a;2 + x + 2) -x + 1 


4.152590 



Table 2-8 shows the comparison of these methods for these functions. 
All the numerical computations have been carried out using MATHE- 
MATICA 8. The numerical results show that the our proposed method 
is efficient. 

Table 2. Errors Occurring in the estimates of the root 
of function fi{x) = sin^(x) — x^ + 1 after third iteration 
by the method described with initial guess xq = 1. 



Methods 



l/i(:^3)| 



Steffensen Method 
Jain Method 
Dehghan I Method 
Dehghan II Method 
Dehghan III Method 
Cordero Method 
MKDF Method 



0.27307e-3 

0.15974e-ll 

0.12208e-3 

0.27058e-7 

0.23473e-9 

0.16292e-9 

0.47200e-25 



NEW EFFICIENT STEFFENSEN TYPE METHOD FOR SOLVING. 



Table 3. Errors Occurring in the estimates of the root 
of function f2{x) = x"^ — e^ — 3x + 2 hj the method 
described with initial guess Xq = 0.7. 



Methods 


\f2{x,) 


Steffensen Method 


0.79151e-6 


Jain Method 


0.99771e-31 


Dehghan I Method 


0.27162e-5 


Dehghan II Method 


0.20573e-22 


Dehghan III Method 


0.70112e-12 


Cordero 


0.89687e-6 


MKDF Method 


0.77899e-31 



Table 4. Errors Occurring in the estimates of the root 
of function f^{x) = cos{x) — a; by the method described 
with initial guess Xq = 1. 



Methods 



l/3(:^3)| 



Steffensen Method 
Jain Method 
Dehghan I Method 
Dehghan II Method 
Dehghan III Method 
Cordero Method 
MKDF Method 



0.82149e-10 
0.10235e-34 
0.17099e-23 
0.61489e-33 
0.19366e-28 
0.14803e-15 
0.10720e-64 



Table 5. Errors Occurring in the estimates of the root 
of function /4(x) = cos(x) — xe^ + x"^ by the method 
described with initial guess xq = 1. 



Methods 



1/4(2:3)1 



Steffensen Method 
Jain Method 
Dehghan I Method 
Dehghan II Method 
Dehghan III Method 
Cordero Method 
MKDF Method 



0.46704e-2 

0.15915e-10 

0.85626e-2 

0.40435e-7 

0.15970e-5 

0.19461e-2 

0.16214e-27 
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Table 6. Errors Occurring in the estimates of the root 
of function /^{x) = e^ — 1.5 — arctan(a;) by the method 
described with initial guess Xq = 1. 



Methods 



l/5(:^3)| 



Steffensen Method 
Jain Method 
Dehghan I Method 
Dehghan II Method 
Dehghan III Method 
CorderoMethod 
MKDF Method 



0.78069e-3 

0.99920e-15 

0.43288e+0 

0.11102e-15 

0.21494e-10 

0.55888e-13 

0.11102e-15 



Table 7. Errors Occurring in the estimates of the root 
of function fQ{x) = 8a; — cos(a;) — 2x'^ by the method 
described with initial guess xq = 1. 



Methods 



1/6(2:3)1 



Steffensen Method 
Jain Method 
Dehghan I Method 
Dehghan II Method 
Dehghan III Method 
Cordero Method 
MKDF Method 



0.77299e+l 

0.37235e-3 

0.56299e-9 

0.29397e-2 

0.52774e-14 

0.77803e-24 

0.38994e-14 



Table 8. Errors Occurring in the estimates of the root 
of function fj{x) = ln(a;^ +x + 2) — x + lhj the method 
described with initial guess Xq = 3.6. 



Methods 



fiixs) 



Steffensen Method 
Jain Method 
Dehghan I Method 
Dehghan II Method 
Dehghan III Method 
Cordero Method 
MKDF Method 



0.51451e-13 
0.11061e-42 
0.15144e-29 
0.48976e-33 
0.19986e-49 
0.43626e-79 
0.39633e-74 
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4. Conclusion 

The problem of locating roots of nonlinear equations occurs fre- 
quently in scientific and engineering works. In this paper, we have 
introduced a new derivative free fourth-order method for solving non- 
linear equations. By theoretical result, we confirm that the new method 
archive fourth order convergence. This method can be used for solving 
nonlinear equations without computing derivatives. The new method is 
also compared in the performance with some well known stefFensen type 
methods. Numerical results indicate that the our proposed method 
gives better performance. 
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